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METHOD AND APPARATUS FOR ACCURACY-AWARE 

ANALYSIS 

Background 

[0001] In general, floating-point notation represents a number in the form 
±0.ddd...d X b^, where Z? is a base, each d is a digit from 0 to b-l, and e is an 
integer exponent. The digits (d) are collectively referred to as a mantissa. The 
number of digits in the mantissa is typically fixed for a particular floating-point 
representation. For example, in the most common current form of floating-point 
representation (Le,, Institute of Electrical and Electronics Engineers (IEEE) 754 
binary standard), the base (b) is 2, the exponent (e) ranges from -1023 to +1024 
(represented using 11 bits), and the mantissa is represented using 53 bits 
(including a bit for the sign of the number). 

[0002] Numbers are represented in floating-point notation in a wide variety of 
computer applications. Specifically, programmers represent numbers in floating- 
point notation when the accuracy of the number is important for the particular 
application. However, a programmer's reliance on the accuracy of a given 
number in floating-point notation may be unfounded. In many cases, simple 
arithmetic operations such as addition, subtraction, multiplication, etc. introduce 
accuracy loss that is often hard to detect, and yet is important to understand. For 
example, the addition of Yz to one million one thousand times, is 0.1000500x10 . 
In a system with seven decimals of precision, due to rounding of the intermediate 
result after each addition of Yi, the result is typically represented as 
0.1000000x10 . Depending on what this particular number, represented in 
floating-point notation, is to be used for, a difference of five hundred may be 
significant. Accordingly, it is desirable to have an understanding of the accuracy 
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of the number represented in the mantissa, when the number is used in a 
computer application. 

[00031 One conventional approach to determine the accuracy of a number 
represented in floating-point notation uses interval arithmetic (lA). Using this 
method, numbers represented using floating-point notation are replaced with 
intervals that bound the result (final or intermediate) when the number cannot be 
expressed exactly by the computer application (i.e., the number of mantissa bits 
required to represented the number exceeds the number of mantissa bits available 
to represent the number, e.g., a number requires nine bits to be exactly 
represented but the mantissa in the given floating-point representation is only 
allocated seven bits). 

[0004] Accordingly, lA provides a means to indicate to the progranuner whether 
the given number, represented in floating-point notation, is exact or inexact. 
However, lA typically provides no information about the cause of the inaccuracy 
or inability of the computer program to accurately represent the result in floating- 
point notation. In addition, lA typically requires low level-access to the rounding 
modes of the processor upon which the computer application is running. 

[0005] A second conventional approach to determine the accuracy of a number 
represented in floating-point notation modifies the level of precision in the 
processor executing the computer application and determines whether the change 
in precision (e.g., from 14 decimals in the mantissa to 11 decimals in the 
mantissa) has an effect on the accuracy of the results. This approach provides a 
means to identify whether a given result is exact or inexact, as well as to 
determine, what level of precision is required to represent the result. However, 
this approach does not typically provide any information as to the cause of the 
inaccuracy of the results. 
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Summary 

[0006] In general, in one aspect, the invention relates to a method for accuracy- 
aware analysis of a program comprising obtaining source code for the program 
comprising a floating point variable, instrumenting the source code to associate 
an accuracy-aware tracking structure with the floating-point variable to obtain 
instrumented source code, compiling to instrumented source code to obtain 
instrumented compiled code, and executing the instrumented compiled code, 
wherein executing the instrumented compiled code comprises using the 
accuracy-aware tracking structure to track an operation on the floating-point 
variable. 

[00071 In general, in one aspect, the invention relates to A system for performing 
accuracy-aware analysis on a program, comprising a source code defining the 
program comprising a floating point variable, a preprocessor configured to 
instrument the source code, wherein instrumented source code comprises an 
accuracy-aware tracking structure associated with the floating-point variable, a 
runtime utility configured to track an operation of the floating-point variable and 
update the accuracy-aware tracking structure associated with the floating-point 
variable, and a runtime environment configured to execute the instrumented 
source code and the runtime utility. 

[0008] In general, in one aspect, the invention relates to A computer system for 
performing accuracy-aware analysis on a program, comprising a processor, a 
memory, a storage device, and software instructions stored in the memory for 
enabling the computer system under control of the processor, to obtain source 
code for the program comprising a floating point variable, instrument the source 
code to associate an accuracy-aware tracking structure with the floating-point 
variable to obtain instrumented source code, compile to instrumented source 
code to obtain instrumented compiled code, and execute the instrumented 
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compiled code, wherein executing the instrumented compiled code comprises 
using the accuracy-aware tracking structure to track an operation on the floating- 
point variable. 

[0009] Other aspects and advantages of the invention will be apparent from the 
following description and the appended claims. 

Brief Description of Drawings 

[0010] Figure 1 shows a flow diagram in accordance with one embodiment of the 
invention. 

[0011] Figure 2 shows a flowchart in accordance with one embodiment of the 
invention. 

[0012] Figure 3 shows a computer system in accordance with one embodiment of 
the invention. 

Detailed Description 

[0013] Specific embodiments of the invention will now be described in detail with 
reference to the accompanying figures. Like elements in the various figures are 
denoted by like reference numerals for consistency. 

[0014] In the following detailed description of embodiments of the invention, 
numerous specific details are set forth in order to provide a more thorough 
understanding of the invention. However, it will be apparent to one of ordinary 
skill in the art that the invention may be practiced without these specific details. 
In other instances, well-known features have not been described in detail to avoid 
obscuring the invention. 

[0015] In general, embodiments of the invention relate to a method and apparatus 
for accuracy-aware analysis. More specifically, embodiments of the invention 
provide a means to track operations performed on number a represented in 
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floating-point notation and identify accuracy loss due to the operations within a 
given computer application. Embodiments of the invention also provide a means 
to modify the source code of the computer application to include functionality to 
track operations performed on number a represented in floating-point notation 
and identify accuracy loss due to the operations within a given computer 
application. 

[0016] Figure 1 shows a flow diagram in accordance with one embodiment of the 

invention. Initially, source code (100), typically in a high-level language {e.g., 
Java™ (Java™ is a trademark of Sun Microsystems, Inc.), Fortran, C, C++, etc.) 
is preprocessed in a preprocessor (102) to produce instrumented source code 
(104). The instrumented source code (104) is functionally identical to the source 
code (100); however, the instrumented source code (104) has been augmented to 
associate an accuracy-aware tracking structure for each number represented in 
floating-point notation (i.e., a floating-point variable). 

[0017] In one embodiment of the invention, the accuracy-aware tracking structure 
records whether the variable is exact or inexact, and if the variable is inexact 
indicates where the uncertainty (e.g., at which bit in the mantissa) of the floating- 
point variable begins. In addition, the accuracy-aware tracking structure may 
record the history of operations that have been performed on the floating-point 
variable. The history of operations may include the number of multiplication 
operations, the number of division operations, the number of square root 
operations, number of addition operations, number of subtraction operations, etc. 
Also, the accuracy-aware tracking structure may also record the number of left- 
digit destructions due to renormalization, number of scaled mantissa digits 
(described below), etc. 

[0018] In one embodiment of the invention, the accuracy-aware tracking structure 
may be implemented as a multi-dimensional array indexed by floating-point 
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variable. Alternatively, the accuracy-aware tracking structure may be 
implemented as a record (such as a struct in C). Further, in one embodiment of 
the invention, each element in an array of floating-point variables may be 
associated with a unique accuracy-aware tracking structure. Those skilled in the 
art will appreciate that the accuracy-aware data structure may be implemented 
using any type of data structure that can store the aforementioned information. 

[0019] Continuing with the discussion of Figure 1, the preprocessor (102), in one 
embodiment of the invention, includes an instrumentation utility (106). The 
instrumentation utility (106) typically includes functionality to associate each 
floating-point variable in the source code (100) with an accuracy-aware tracking 
structure. The instrumentation utility (106) typically uses a lexical analyzer 
(108) to parse the source code (100) to obtain tokens, and a parse generator (110) 
to identify a sequence of tokens. Upon identification of the sequence of tokens, 
the instrumentation utility (106) inserts the appropriate code into the source code 
(100). 

[0020] In general, the code inserted into the source code (100) functions as a hook 
to record the occurrence of certain events within the source code (100). More 
specifically, in one embodiment of the invention, the code inserted within the 
source code (100) may include functionality to call into a runtime logging utility 
(118) (described below), functionality to associate an accuracy-aware tracking 
structure with a floating-point variable, functionality to update the accuracy- 
aware tracking structure associated with the floating-point variable after an 
operation has been performed on the floating point variable, etc. Those skilled in 
the art will appreciate that while the instrumentation utility (106) has been 
described as inserting code into the source code (100), the instrumentation utility 
(106) may also include functionality to modify the source code (100), as opposed 
to strictly inserting code. 
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[0021] Once the instrumented source code (104) has been generated, the 
instrumented source code (104) is compiled, via a compiler (112), to generate 
instrumented compiled code (116). The compiler (112) typically links the 
instrumented source code (104) with the necessary libraries to support the 
functionality of the code inserted into the source code (100) to generate the 
instrumented source code (104). For example, in some embodiments of the 
invention, the instrumented source code (104) may include calls into the runtime 
logging utility (118) (described below). Thus, when the instrumented source 
code (104) is compiled, the compiler (112) links the libraries associated with the 
runtime logging utility (118) with the instrumented compiled code (116). 

[0022] The instrumented compiled code (116) is subsequently executed in a 
runtime environment (114), for example in a Java™ Virtual Machine. The 
runtime environment (114) typically includes functionality to execute the 
instrumented compiled code (116) as well as the functionality to implement 
additional utilities linked to the instrumented compiled code (116). For example, 
the runtime environment (114) may concurrently execute a runtime logging 
utility (118). In one embodiment of the invention, the runtime logging utility 
(118) includes functionality to maintain the accuracy-aware tracking structure 
and functionality to update the accuracy-aware tracking structure associated with 
a given floating-point variable after an operation has been performed on the 
floating-point variable. 

[0023] Further, in one embodiment of the invention, the runtime logging utility 
(118) also includes functionality to perform the necessary operations to update 
the accuracy-aware data structures (e.g., determine the scaled mantissa 
(described below)). In addition, the runtime logging utility (118) may also 
include functionality to generate an accuracy report (122) for one or more 
floating-point variable within the source code (100). Note that once the 
instrumented compiled code (116) has finished executing, the runtime 
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environment (114) may also display the results (120) (i.e., values) of the floating- 
point variables. 

[0024] The following discussion details an embodiment of the sequence of events 
that occur when the instrumented compiled code (116) is executed in the runtime 
environment (114). Figure 2 shows a flowchart in accordance with one 
embodiment of the invention. Specifically, Figure 2 shows a method for 
accuracy-aware analysis of a computer program. Initially, an operation of on a 
floating-point variable is performed (Step 200). A hook, provided by the inserted 
code, is subsequently triggered that calls into the runtime logging utility (Step 
202). The runtime utility subsequently updates the operation counter, within the 
associated accuracy-aware tracking structure, corresponding to the operation 
performed on the floating-point variable. For example, if an addition operation 
was performed on a floating-point variable, then the addition variable would be 
incremented to indicate the occurrence of the addition operation. Altematively, 
depending on whether the operation performed on the floating-point variable 
included renormalization, a variable tracking additions with renormalization or a 
variable tracking additions without renormalization may be incremented. 

[0025] Continuing with Figure 2, the runtime logging utility subsequently 

determines a scaled mantissa for each floating-point variable present in the 
operation (i.e., the operands), as well as the resulting floating point variable (Step 
206). The scaled mantissa representation for a floating-point variable 
corresponds to the integer represented by the first non-zero digit in the mantissa 
to the last non-zero-digit in the mantissa. Thus, for example, the scaled mantissa 
for the number -0.1 1 1001001 10010000000 x 2^^ is -0.1 1 1001001 1001. 

[0026] Continuing with Figure 2, the scaled mantissa for the result is subsequently 

used to determine whether the result is exact or inexact (Step 208). More 
specifically, if the floating-point variable operands are assumed to be exact, then 
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the result of the operation may be exact or inexact. In one embodiment of the 
invention, the result is said to be inexact when the correct result, expressed as a 
scaled mantissa, requires more bits to represent than there are in the mantissa of 
the format for the result. For example, if a binary format allows a four-bit 
mantissa, then (0.1110 x 2^) (0.1010 x 2^) is exactly expressible as 0.100011 x 
2^; the scaled mantissa for the result is 110011. Accordingly, because 110011 
requires six bits to express but only four bits of mantissa are available, the result 
is rounded to 0.1001 x 2^, and marked as inexact. 

[0027] Inexact results are similarly produced when addition, subtraction, division, 
or any other mathematical operation on quantities represented as floating-point 
values produce a result that does not fit in the fixed mantissa size. Note that 
every high-level language has an operation for finding the remainder of a 
division, so in the case of a division operation, if the remainder is zero the 
resulting quotient is exact, otherwise the result is flagged as inexact. 

[0028] In one embodiment of the invention, the inexact condition is associated 
with the floating-point variable and stored in the associated accuracy-aware 
tracking structure. In such cases, the inexact condition may be represented by a 
one-bit flag. Alternatively, the inexact condition may be associated with the 
floating-point variable and stored in a temporary memory store that is separate 
from the accuracy-aware tracking structure and that persists throughout the 
duration of the execution of the instrumented compiled code. 

[0029] As an altemative to determining the scaled mantissa, embodiments of the 
invention using the IEEE 754 standards may use a condition flag incorporated 
within the standard. Note that the condition flag specified in the IEEE 754 
standard requires the runtime logging utility (or a related process) to have access 
to the "sticky bit" which represents the flag, and functionality to maintain 
knowledge of the inexact condition after the result has been computed. 
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[0030] Returning to the discussion of Figure 2, if the resuh is determined to be 
exact (Step 206), then the number of half units in last place (HULP) associated 
with the result is not incremented (Step 210). In one embodiment of the 
invention, the term HULP is used to record the error caused by rounding the last 
bit(s) in the mantissa, error introduced by renormalization (typically due to left 
digit destruction), etc. Thus, the number of HULP associated with a given 
operation corresponds to the base of the floating-point representation raised to 
the power of the number of bits causing the error. 

[0031] For example, consider the subtraction of a first floating-point variable from 
a second floating point- variable: 

0.3141593 X 10^ 
0.3140000 X 10^ 
= 0.0001593x 10^ 

The result, 0.0001593 x lOMs subsequently renormalized to 0.1593000x lO'l As 
shown above, the renormalization of the result destroys the left three zeros after 
the decimal place, and introduces three zeros to the right of the last non-zero 
integer in the mantissa. Thus, the number of decimals of error introduced into the 
results is three, and accordingly, the number of HULP recorded for this operation 
is 10^ or one thousand. The HULP variable in the accuracy-aware tracking 
structure is subsequently incremented by one thousand. Thus, the invention 
informs the programmer or user that the result is only accurate to the fourth 
decimal place. Further, the invention provides the history of operations that 
caused the inaccuracy. Note that if the number is inexact and no renormalization 
of the result has occurred, then the number of HULP associated with the operation 
is the base raised to the power of zero, Le,, the HULP is one. 

[0032] Those skilled in the art will appreciate that while HULP are described as a 
means to quantify the error associated with a given number represented in 
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floating-point notation, the invention is not limited to quantifying error using 
HULP. For example, as an altemative to using HULP to quantify error, the 
invention may quantify error by using intervals. Consider the following 
example: If one were to express the exact number 33 as the inexact number 
0.3333x10 , the number of HULP would be 1 (assuming a four-decimal 
mantissa). However, using intervals, the error may be quantified as 0.3333x10^ 
+ (0.00003333x10^ to 0.00003334x10^), where the interval (0.00003333x10^ to 
0.00003334x10^) is used to quantify the error associated with number 
represented in floating-point notation. 

[0033] In one embodiment of the invention, the numbers that make up the interval 
are stored as two separate variables, Le,, the number corresponding to the upper 
limit of the interval and the number corresponding to the lower limit of the 
interval. Further, the numbers that make up the interval may also be represented 
as floating-point variables. Similar to HULP, intervals may be tracked and 
updated after every operation on a floating-point number. Additional 
mechanisms for tracking and quantifying error may also be implemented within 
the scope of the invention. 

[0034] Continuing with Figure 2, those skilled in the art will appreciate that step 

(204) may not need to be performed if at least one of the operands used to 
compute the result was inexact. In either case, whether two exact floating-point 
variables generated an inexact floating-point variable or whether one or more of 
the operands was inexact, the number of HULP for the particular operation is 
determined (Step 212). 

[0035] If the number of mantissa bits required to represent the result exceeds the 
number of bits allocated to represent the mantissa and no renormalization of the 
result occurs, then the number of HULP for the operation is one. However, if 
renormalization of the result occurs, then the number of HULP equals the base 
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raised to the power of the number of zeros added to the mantissa after the last 
non-zero integer in the result. The HULP counter associated with the floating- 
point variable containing the result is subsequently updated using the calculated 
HULP for the renormalization (Step 214). 

[0036] In one embodiment of the invention, the HULP is represented as an integer. 
Specifically, the HULP may be represented by an integer of the same word size 
as the floating-point variable. For example, a 32-bit floating-point variable with 
a 24-bit mantissa may lose 2 HULP before the entire mantissa becomes invalid, 
and 2 fits into a 32-bit integer. Integers of the aforementioned size are typically 
available in processor designs as a native type, or they are readily maintained 
with software that extends the precision of the largest integer native type of the 
system. 

[0037] At this stage, after the HULP for the operation has been determined and the 
appropriate HULP counter has been updated, the runtime logging utility may 
return control back to the instrumented compiled code to continue executing 
within the runtime environment. However, in some embodiments of the 
invention, prior to returning to executing the instrumented compiled code, the 
runtime logging utility (or a related process) determines whether a given 
accuracy threshold has been exceeded with respect to the result (Step 214). If the 
accuracy threshold is not exceeded, then the instrumented compiled program 
continues executing (Step 216). Alternatively, if the accuracy threshold is 
exceeded, then the runtime logging utility (or a related process) performs the 
appropriate action as specified by, for example, the user (Step 218). In one 
embodiment of the invention, the appropriate action to perform is to halt the 
execution of the instrumented compiled code. 

[0038] In one embodiment of the invention, the determination of whether the 
accuracy threshold has been exceeded uses relative error. Altematively, the 
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determination of whether the accuracy threshold has been exceeded uses absolute 
error. Furthermore, the determination of whether the accuracy threshold has 
been exceeded may be that the original program performs a comparison test, 
(e.g., "if a<b then...") on variables which are so approximate that their ranges 
of possible values overlap. For example, if a=0.42 ± 0.2, and Z?=0.43±0.2, then 
the truth of a<b cannot be determined and may be further used to halt execution. 

[0039] Consider the following example, ax bx c = 0^ where the roots are 
expressible as x = {-b ± Qp" - Aacf^^) I {2a). The vast majority of today's 
programmers would program this closed-form expression as written above, using 
and case of the to give the two roots (or one degenerate root ifb- 
4ac is exactly zero.) This approach is hazardous to accuracy, because (b^ - 
Aacf^^ might be quite close to b, and then -b + {b^ - Aacf^^ would cause 
destructive cancellation of digits in the result. 

[0040] Continuing with the above example, suppose a = 3, 6 = 100, and c = 2, that 
the values of a, Z?, and c are exact, and that the mantissa includes 7 digits. A 
program to solve for the root may then perform the following operations: 



tl=b*b 




0.1000000 X 10^ 


t2=4*a 


4a 


0.1200000 X 10^ 


t2=t2*c 


4ac 


0.2400000x10^ 


t2=tl-t2 


b^-Aac 


0.9976000x10* 


ift2 0, then print 


Exit if solution is 




"degenerate or 


degenerate or involves 




non-real answer" 


imaginary numbers 




stop 






end if 






t2=sqrt(t2) 




0.9987993 x 10^ 


rl=-b+t2 


-b^ib^-Aacf^ 


0.1200700 X 10° 


rl=rl/2 


(-b+(b^-4acfy2 


0.6003500 X 10"' 


rl=rl/a 


First root 


0.2001 167 X lO*' 


r2=-b-t2 


-b-ib^-Aacf^ 


-0.1998799x10^ 
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r2=r2/2 -0.9993995 x 10^ 

r2=r2/a Second root -0.3331332 x 10^ 

print rl, r2 
end 

[0041] The result of executing the above algorithm provides two roots, namely 
0.2001167 and -33.31332, to seven decimals. At this stage, there is no 
indication of the accuracy of the two roots. Using the invention, as described 
above, the following accuracy report may be generated. The following table 
shows columns for the number of HULP, the number of Scaled Mantissa digits 
("S.M."), the computed value, the number of addition/subtraction operations 
where the result is not renormalized "±", the number of left digit destructions that 
require renormalization of the result "Idd", the number of multiplication "*", 
division "/", and square root " "operations. The number of S.M. digits is not the 
same as the number of significant digits, so the number 100 may have exactly 
three significant digits, but only one S.M. digit in a decimal floating-point 
representation, the scaled mantissa is the mantissa stripped of leading zeros and 
trailing zeros, without regard for their significance. 





HULP S.M. 


Value 


± Idd 


* 


/ 




tl=b*b 


0 


1 


0.1000000x10^ 


0 


0 


1 


0 


0 


t2=4*a 


0 


2 


0.1200000 X 10^ 


0 


0 


1 


0 


0 


t2=t2*c 


0 


2 


0.2400000 X 10^ 


0 


0 


2 


0 


0 


t2=tl-t2 


0 


4 


0.9976000 X 10^ 


1 


0 


3 


0 


0 


ift2 0, then print 




Test is safe; 














"degenerate or 




no ambiguity 














non-real answer'' 


















stop 


















end if 


















t2=sqrt(t2) 


1 


7 


0.9987993 x 10^ 


1 


0 


3 


0 


1 


rl--b+t2 


100 


5 


0.1200700 X 10^ 


1 


1 


3 


0 


1 


rl=rl/2 


100 


5 


0.6003500 X 10"' 


1 


1 


3 


1 


1 


rl=rl/a 


101 


5 


0.2001167 X 10° 


1 


1 


3 


2 


1 


r2=^b-t2 


2 


7 


-0.1998799 X 10^ 


2 


0 


3 


0 


1 


r2=r2/2 


3 


7 


-0.9993995 x 10^ 


2 


0 


3 


1 


1 
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r2=r2/a 4 7 -0.3331332x 10^ 2 0 3 2 1 

print rl, r2 

end 

[0042] Upon inspection of the above report, it is clear that the first root (rl) is only 
accurate to 5 decimal places while the second root (r2) is accurate to 7 decimal 
places. 

rl=rl/a 101 5 0.2001167x 10^ 113 2 1 

r2=r2/a 4 7 -0.3331332x 10^ 2 0 3 2 1 

[0043] Referring to the above report, the annotations indicate that the rl 
computation involved a destructive subtraction (Idd = 1), whereas the r2 
computation did not (Idd = 0). Specifically, the rl computation involved adding 
-100 to +99.87993, after which the result was renormalized by shifting the result 
left two decimal places. That renormalization incremented the Idd counter by 1 
and the corresponding HULP counter by 100. 

[0044] Using the above annotations, the programmer may reformulate the 
algorithm to avoid the destructive subtraction. For example, the following 
modification to the above algorithm may restore the accuracy of rl: Compute r2 
first, then compute rl=c/(a*r2). Note that in one embodiment of the invention 
the annotations indicate where reformulations are needed. In addition, the 
annotations may also suggest to the programmer to use interval arithmetic for 
containing the error or to use higher-precision variables for part of the 
computation. Those skilled in the art will appreciate that the above example was 
intended to illustrate the operation of the invention, but was not intended to limit 
the scope of the invention. 

[0045] In one embodiment of the invention, the invention may be used to 
determine locations within the source code where using a lower level of precision 
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would still produce accurate results. For example, the invention may be used to 
identify locations in the code where 32-bit floating-point arithmetic could be 
used to store temporary values without affecting the accuracy of the final result 
that in other places relies on 64-bit floating-point arithmetic. This type of 
optimization may be used to effectively reduce both the storage and bandwidth 
requirements for a given algorithm at no penalty to the user. Further, if the 
invention provides enough confidence that the lower precision is safe for the 
range of the possible data inputs in an a given computer application, then a 
computer application may be modified to calculate results at a faster rate while 
using less memory resources to store the floating-point variables without 
reducing the accuracy of the numerical results. 

[0046] Those skilled in the art will appreciate that the invention may be 
implemented on any high level language, e.g., Java™, Fortran, C, C++, etc., that 
supports the use of floating-point calculations. Further, the invention typically 
does not require any extension to the high level language, and typically does not 
violate any security or portability specifications. Further, those skilled in the art 
will appreciate that while the invention has been described using a given 
floating-point representation, the invention may be implemented using any 
floating-point representation. 

[0047] The invention may be implemented on virtually any type of computer 
regardless of the platform being used. For example, as shown in Figure 3, a 
networked computer system (300) includes a processor (302), associated memory 
(304), a storage device (306), and numerous other elements and functionalities 
typical of today's computers (not shown). The networked computer (300) may 
also include input means, such as a keyboard (308) and a mouse (310), and 
output means, such as a monitor (312). The networked computer system (300) is 
connected to a local area network (LAN) or a wide area network {e.g., the 
Intemet) (not shown) via a network interface connection (not shown). Those 

16 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 03226.367001 ; SUN040407 

skilled in the art will appreciate that these input and output means may take other 
forms. Further, those skilled in the art will appreciate that one or more elements 
of the aforementioned computer (300) may be located at a remote location and 
connected to the other elements over a network. 

[0048] In one or more embodiments of the invention, a means to readily automate 
a given program to track operations performed on numbers represented in 
floating-point notation is provided. Also, an indication as to the possible sources 
of inaccuracy is provided. Further, embodiments of the invention provide a 
means to pinpoint an operation where accuracy drops precipitously and to notify 
the programmer of such as occurrence. Embodiments of the invention also 
provide insight into the type of accuracy loss resulting from small losses of 
accuracy that accumulate with repetition to affect the results. By tracking 
operations performed on numbers represented in floating point notation, 
embodiments of the invention provide a means to increase programmer 
productivity and decrease the amount of time required to currently find such 
errors. 

[0049] While the invention has been described with respect to a limited number of 
embodiments, those skilled in the art, having benefit of this disclosure, will 
appreciate that other embodiments can be devised which do not depart from the 
scope of the invention as disclosed herein. Accordingly, the scope of the 
invention should be limited only by the attached claims. 
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